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In this work we consider Randall-Sundrum brane-world type scenarios, in which the spacetime 
is described by a five-dimensional manifold with matter fields confined in a domain wall or three- 
brane. We present the results of a systematic analysis, using dynamical systems techniques, of 
the qualitative behaviour of the Friedmann-Lemaltre-Robertson-Walker and the Bianchi I and V 
cosmological models in these scenarios. We construct the state spaces for these models and discuss 
how their structure changes with respect to the general-relativistic case, in particular, what new 
critical points appear and their nature, the occurrence of bifurcations and the dynamics of anisotropy. 
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I. INTRODUCTION 

String and membrane theories are promising candi- 
dates for a unified theory of all forces and particles in 
Nature. A consistent construction of a quantum string 
theory is only possible in more than four spacetime di- 
mensions. Then, to make a direct connection of these 
theories with our familiar non-compact four-dimensional 
spacetime we are compelled to compactify the extra spa- 
tial dimensions to a finite size or, alternatively, find a 
mechanism to localize matter fields and gravity in a lower 
dimensional submanifold. 

Recently, Randall and Sundrum have shown that for 
non-factorizable geometries in five dimensions there ex- 
ists a single massless bound state confined in a domain 
wall or three-brane This bound state is the zero 
mode of the Kaluza-Klein dimensional reduction and cor- 
responds to the four-dimensional graviton. The picture 
of this scenario is a five-dimensional Anti-de Sitter space 
{bulk) with an embedded three-brane where matter fields 
are confined and Newtonian gravity is effectively repro- 
duced at large-scale distances. Earlier work on Kaluza- 
Klein dimensional reduction and matter localization in a 
four-dimensional manifold of a higher-dimensional non- 
compact spacetime can be found in 0. 

The Randall-Sundrum model was inspired by string 
theory. In the context of dimensional reduction of eleven- 
dimensional supergravity, Hofava and Witten showed 
that the ten-dimensional Eg x Eg, heterotic string is 
connected with an eleven-dimensional theory compact- 
ified on the orbifold R 10 x S 1 /Z 2 f§. Moreover, they 
concluded that the coupling constants of gauge fields 
in the ten-dimensional boundary are related with the 
eleven-dimensional gravitational constant Q. The pic- 
ture coming out of this model is that of two separated 
ten-dimensional manifolds. Gauge fields are confined in 
these boundary manifolds whereas gravity can propa- 
gate in the higher dimensional spacetime. As a conse- 
quence, these two separated worlds can only communi- 



cate through gravitational interactions. The cosmolog- 
ical implications of the Hofava- Witten theory have al- 
ready been extensively analyzed g. 

The original motivation for the Randall-Sundrum 
model was the solution of the hierarchy problem in a 
slightly different set up || . In this case one has two par- 
allel branes with opposite tensions embedded in a five- 
dimensional spacetime with negative cosmological con- 
stant. Actually, the fifth dimension is compactified in 
the orbifold S 1 /Z2 and the two branes are located at the 
singular boundary points. Due to an exponential fac- 
tor in the metric tensor, the particles living in the neg- 
ative tension brane acquire effectively a huge physical 
mass parameter compared to the fundamental scale with 
a moderate fine tuning of the size of the extra dimen- 
sion. Unfortunately, as pointed out in 0, the cosmology 
in this brane is rather unsatisfactory because the energy 
density of matter present in the brane must be negative, 
which violates the weak energy condition. The argument 
is based on the observation made by Binetruy et al. 
(see also JlO(| ) that the effective Friedmann equation for 
the Hubble parameter for a five-dimensional spacetime 
with energy density localized in a infinitely thin domain 
wall is modified with respect to the general relativistic 
case. Other attempts to solve the hierarchy problem 
in the context of extra dimension have been examined 



in 111). 

Nevertheless, the model with a non-compact fifth di- 
mension and only one brane is consistent with present 
gravity experiments. In general, scenarios with extra 
dimensions predict corrections to the Newtonian poten- 
tial at short distances and important deviations from the 
standard evolution of the universe at early times. Then, 
current day cosmological observations, such as the age of 
the universe or the abundances of light elements, can- 
not be used to constraint these models. In contrast, 
the search for deviations of Newton's law is their fun- 
damental observational probe jl2) . The fact that Newto- 
nian gravity has been tested quite accurately up to 1 mm 
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(M* ~ 10 _16 TeV) limits the value of the fundamental 
scale associated with the hve-dimensional gravitational 



10 5 TeV 



coupling constant: M (5) P > NlJ M, 
p|,^3| (Mp is the Planck mass). Future experiments will 
further constraint this naive estimate 1 14 lj] 

The purpose of the present work is to study the cos- 
mological evolution of these brane-world scenarios. We 
are going to follow the geometric formulation and gen- 
eralization of the Randall- Sundrum scenario introduced 
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The Einstein equations in the bulk can be 
written in the following formjj 



G 



(5) 
AB 



with 



T 



AB = S (x) [-A.9AB +T Al 



(1) 



(2) 



In these expressions Kr^) is the hve-dimensional gravita- 
tional coupling constant; g AB , G^ B and A( 5 ) are the met- 
ric, Einstein tensor and the negative cosmological con- 
stant of the bulk spacetime, respectively; Tab is the mat- 
ter energy-momentum tensor; the spacelike hypersurface 
x 4 = x = gives the brane world and §ab is its induced 
metric; finally, A is the tension of the brane, which must 
be assumed to be positive in order to recover conven- 
tional gravity on the brane. Using the Gauss-Codacci 
equations relatingthe four- and five-dimensional space- 
times, equations ([j}]!) l ea d to the following modification 
of the Einstein's equations of General Relativity on the 
brane (u||l7j 



G a b 



-Ag Q b + n 2 T ab + K( 5 jS a b 



E 



(5) 



(3) 



where g ao is the four-dimensional metric on the brane 
and G a b its Einstein tensor. The four-dimensional gravi- 
tational constant k and the cosmological constant A are 
given in terms of the fundamental constants in the bulk 
by 



= — Xk 



A = 



|A(5)| 



(5) 



respectively; where A c is a critical brane tensior 
by 



given 



* Upper-case Latin letters denote coordinate indices in the 
bulk spacetime (A, B, . . . = 0, . . . , 4) whereas lower-case 
Latin letters denote coordinate indices in the four-dimensional 
spacetime where matter is confined (a, b, . . . = 0, . . . , 3). We 
will use physical units in which c = 1. 

^The particular Randall- Sundrum solution corresponds to 
the case when the tension of the brane A equals the critical 



brane tension A c 



and T ab = E, 



(5) 



: 0. 



K = 6 



|A (5 



"(5) 



(4) 



S ao are corrections quadratic in the matter variables (due 
to the form of the Gauss-Codacci equations) and given 
by 



S a b 



— TT 

12 1 ±a 



\T„ C T, 



be 



±g ab [3T cd T cd - T 2 ] , (5) 



(5) 

where T = T a a - And finally, E ab are corrections coming 

(5) 

from the extra dimension, more precisely, E ab are the 
components of the electric part of the Weyl tensor of the 
bulk, C A 5 BCD , with respect to the normal, n A (n A riA = 
1), to the hypersurface x = where matter is confined, 
that is, 



K (5) _ r (5) C D 

^AB — ^ ACBD n 



Moreover, it is worth to note that the twice contracted 
second Bianchi identities in the bulk, G AB — 0, im- 
ply 



V a T a b 



0, 



(6) 



where we have taken jv, x a } to be Gaussian normal co- 
ordinates (see, e.g., [f9[ ) adapted to the hypersurface 
X = 0. Therefore, we can say that the Einstein equa- 
tions in the bulk (gj) imply the conservation of the energy- 
momentum tensor in the brane world. 

In this paper we will deal with generalized Randall- 
Sundrum scenarios in which the effects of the extra- 
dimension come from the term quadratic in the energy- 
momentum tensor, i.e. S a b (||)- Thus, we are assuming 



E 



(5) | 

ab\x=o 







E 



(5) 



0. 



(7) 



This includes conformally-flat bulks (Cabcd = 0), 
and in particular, the five-dimensional Anti-de Sitter 
spacetime, the bulk considered in the original Randall- 
Sundrum scenario. The extension of this work to general 
bulks will be presented in a future paper [pof . 

For the scenarios just outlined we have constructed 
and studied the state space of the Friedmann-Lemaitrc- 
Robertson- Walker (FLRW) and the Bianchi I and V cos- 
mological models. Then, we have discussed systemati- 
cally how the extra dimension changes the dynamics with 
respect to the general-relativistic case. In particular, we 
find a new critical point representing the dynamics at 
very high energies, in the early universe (near the Big- 
Bang) and also near the Big-Crunch in the case of rec- 
ollapsing models. We also find new bifurcations in the 
state space as the equation of state of matter changes 
(we will assume a perfect-fluid energy-momentum con- 
tent), which are characterized by the occurrence of an 
infinite number of non-general-relativistic critical points. 
Finally, the Bianchi I and V models will provide informa- 
tion regarding the dynamics of anisotropy in the branc- 
world scenario. 
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The paper is organized as follows. In section |l| we 
will study the dynamics of the FLRW models in the 
brane-world scenarios, introducing the notation and some 
tools used in the analysis of dynamical systems (see, 
e.g., [^l],|2)). In section H, we will study the dynam- 
ics in homogeneous but anisotropic cosmological models. 
In particular, the dynamics of the orthogonal Bianchi I 
and V cosmological models, which contain the flat and 
the negatively curved FLRW models, respectively. We 



will finish with some concluding remarks in section IV 



II. 



DYNAMICS OF THE FLRW MODELS IN THE 
BRANE-WORLD SCENARIO 



In this section, we start by assuming that the brane- 
world is described by a FLRW metric. The FLRW space- 
times are the standard cosmological models. As is well- 
known p3p^ ], they are motivated by the so-called Cos- 
mological Principle in the sense that they are homoge- 
neous and isotropic cosmological models (they have a 
six-dimensional group of motions). Then, the line ele- 
ment in the brane-world (x = 0) will be given by 



ds 2 



where 



-dr 



a 2 {t) [dr 2 + E 2 k (r){d6 2 + sin 2 9dtp 2 )] 



sin r for k = 1 , 
E fc (r) = { r for k = 0, 
sinh r for k = — 1 , 

and a(t) is the scale factor. 

Here, we will study the dynamics of the FLRW models 
considering a bulk spacetime satisfying the condition (0), 
which includes the five-dimensional Anti-de Sitter space- 
time. On the other hand, we will assume that the matter 
content is equivalent to that of a perfect fluid and there- 
fore, the energy-momentum tensor will have the following 
form 

Tab = {p + p)u a u b + pg ab , 

where u, p and p are the unit fluid velocity of matter 
{u a u a — —1), the energy density and the pressure of the 
matter fluid respectively. We will also assume a linear 
barotropic equation of state for the fluid, that is, 



P= (7- !)/>• 



(8) 



The weak energy condition (see, e.g., |23|1 ) imposes the 
restriction p > 0, and from causality requirements, the 
speed of sound [c s = (dp/dp) 1 / 2 } must be less than the 
speed of light, we have that 76 [0, 21 . Then, taking into 
account the form of the equations (||,|5|) and (]?]) , it turns 
out that the fluid velocity u is aligned with the velocity of 
the preferred observers in the FLRW spacetimes (except- 
ing in the case G ao oc g ab , where there are no preferred 
observers), those that observe the matter distribution to 



be homogeneous and isotropic. Then, we can write u as 
follows 



u 



d_ 
dt 



u = dt . 



Finally, taking into account recent observations |25| , |26| , 
we will consider only the case of a positive cosmological 
constant, i.e. A > 0. Then, introducing the Hubble 
function H(t) 

a dt a 

the dynamics of the FLRW models imposed by the 
modified Einstein field equations (0) and the energy- 
momentum conservation equation (|6|) is governed by the 
following set of ordinary differential equations 



H 



-H 2 



37-2 



K 2 p 



1 



37- 1 P 

37- 2 A 



p = SjHp, 



+ ^A, (9) 



(10) 



H 



1 



-K'p 1 



6 



-A. 



(11) 



where 3 i? denotes the scalar curvature of the hypersur- 
faces orthogonal to the fluid velocity, the {t = constant} 
hypersurfaces, which is given by 3 R = 6ka~ 2 (t). Equa- 
tion (0) is the modified Raychaudhuri equation, equa- 
tion (|0|) comes from the energy-momentum tensor con- 
servation equation, and finally, equation ( pd| ) is the mod- 
ified Friedmann equation. As is well-known, (|^) is a con- 
sequence of (|l(]) and ( p"l| ) , and the dynamics is completely 
described by the functions (H, p) and the parameters k, 
K, 7, A and A. 

In order to study the dynamics of these models we 
will closely follow the analysis carried out by Goliath & 
Ellis p7| for general relativistic FLRW models with a 
cosmological constant. To that end, and in order to get 
compactified state spaces, it is convenient to consider two 
differentiated cases: (i) 3 R < (k — or k — — 1) and 
(ii) 3 R > (fc = 1). 

In the case (i), let us introduce the following set of 
dimcnsionless variables 



Sir. 



SI, 



2 

up 
ZIP 
A 
MP 



Slk = 



[i R _ 
6H 2 ~ 

K 2 p 2 



k 



1 



6A H 2 



(12) 
(13) 



where J7 P is the ordinary density parameter and fife, Q\ 
and Q\ are the fractional contributions of the curvature, 
cosmological constant and brane tension, respectively, to 
the universe expansion ( fil| ) . Therefore, all of them have 
a clear physical meaning. As we can see, they are non- 
negative and singular when H = 0. Furthermore, the 
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Friedmann equation (11), which now takes the following 
simple form 



n n + n k + riA + fi> 



(14) 



the interval [0, 
coordinates fi 
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implies that they must belong to 
and hence, the state space with 
(fip, fife, fiA, is compact. 

In order to find the dynamical equations for these vari- 
ables we will introduce the following dimensionless time 
derivative 



1 d 
W\di' 



(15) 



where \H\ is the absolute value of H. Then, we have 



H' = -e(l + q)H, 



(16) 



where e is the sign of H [e = sgn(H)]. As is clear, for 
e = 1 the model will be in expansion, and for e = — 1 
it will be in contraction. Moreover, q is the deceleration 
parameter, which is defined by 



1 a 



3 7 - 2 



fip - 



(37-i)n> 



Then, the dynamical system for our dimensionless vari- 
ables (|l2],[l3|) can be written in the following form 



fip = e[2(l + <?)-3 7 ] n p , 

fi' fc = 2egfi fe , 

fi A = 2e(l + g)fi A , 

n' x = 2e(l + q-3 1 )n x . 



(17) 
(18) 
(19) 
(20) 



It is important to note that equation (|l^) is not coupled 
to the system of equations ( p"7fp0| ) , and therefore we can 
ignore it for the dynamical analysis. To begin with, we 
have to find the critical points of this dynamical system, 
which can be written in vector form as follows 



n' = /(n), 



where / can be extracted from (JT7-20). The critical 
points, fi*, which are the points at which the system 
will stay if initially it was there (see, e.g. p2[| )) are given 
by the condition 

/(n*) = o. 

Their dynamical character is determined by the eigenval- 
ues of the matrix 



on 



fi=rr 



If the real part of the eigenvalues of a critical point is not 
zero, the point is said to be hyperbolic. In this case, the 
dynamical character of the critical point is determined 
by the sign of the real part of the eigenvalues: If all 



of them are positive, the point is said to be a repeller, 
because arbitrarily small deviations from this point will 
move the system away from this state. If all of them 
are negative the point is called an attractor because if we 
move the system slightly from this point in an arbitrary 
way, it will return to it. Otherwise, we say the critical 
point is a saddle point. The dynamical system ( |l7}]2C| ) 
has four hyperbolic critical points corresponding to: the 
flat FLRW models (F), k = A = A" 1 = and a(t) = 
1 2/(3 7 ). the Milne universe (M), p = A = 0,fe = -l and 
a(t) = t; the de Sitter model (dS), k = p = and a(t) = 
exp( 1k /A/3i); and a non-general-relativistic model (m) 
first discussed by Binetruy, Deffayet and Langlois 
brane-world scenario without brane tension (see ^, 
more details). Their coordinates in the state space, i.e., 
fi = (dp, fife, Oa, £7a), and their eigenvalues are given in 
the following table p8| 




Model 


Coordinates 


Eigenvalues 


F e 


(1,0,0,0) 


e(3 7 - 2, 3 7 - 2, 3 7 , -3 7 ) 


Me 


(0,1,0,0) 


e(-(3 7 - 2), 0,2, -2(37- 1)) 


dS e 


(0,0,1,0) 


-e(3 7 ,2,2,6 7 ) 


m c 


(0,0,0,1) 


6(37,2(37-1), 67, 2(37-1)) 



The dynamical character of these points is given in a 
table below. 

Now, let us consider the situation of the case (ii), in 
which 3 R is positive. As we have already mentioned, 
in this case the state space defined by the variables 
11 = (fip, fife, fiA, ^a) is no longer compact (because 
now fife < 0). However, we can introduce another set 
of variables, analogous to the ones introduced previ- 
ously ( |l2| , [13| ), describing a compact state space. Firstly, 
instead of using the Hubble function H we will use the 
following quantity 



D = 



(21) 



and from it, let us define the following dimensionless vari- 
ables 



fi. 



H 
= D 
A 
3~D2 



fi =^ 
p ~ 3D? 



(22) 



fiA = 



1 K 2 P 2 

6A D 2 



From these definitions we see that now the case H = 
is included. Moreover, the Friedmann equation takes the 
following form 



fi A + fiA = 1 ■ 



(23) 



which, together with the fact that — 1 < Q < 1 [see 
Equation (p2|)], implies that the state space defined by 
the new variables is indeed compact. Using the following 
new time derivative 



)_d_ 

Ddi ' 



(24) 
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the system of evolution equations for the variables D, Q, 
fl p , ft a, and fl\ is given by 

D' = -{l + qQ 2 )QD, 

Q' = -qQ 2 (l-Q 2 ), (25) 

Sl' p = [2(l + qQ 2 )-3j]QCl p: (26) 

Sl A = 2(1 + qQ 2 )QSl\ , (27) 

Sl' x = 2 [1 + qQ 2 - 3 7 ] QSl x , (28) 

where the deceleration parameter is now given by 
1 + qQ 2 = ^(Sl p + 2Sl x ), 

The evolution equation for D is not coupled to the 
rest, so we will not consider it for the dynamical study. 
Thus, we study the dynamical system for the variables 
SI = (Q,Sl p ,Sl\,Sl x ), determined by the equations ( p5[ - 
p8| ). The complete set of critical points, their coordinates 
in the state space, i.e. SI* , and their corresponding eigen- 
values are given in the following table p8| 



Model 



Coordinates 



Eigenvalues 



dS e 
E 

111,: 



(e, 1,0,0) 
(e, 0,1,0) 

(e, 0,0,1) 



£(37-2,37,37,-37) 
-6(2,87,0,67) 

(0, v&o.-v^) 

2e(3 7 - 1,^,37, 3 7 ) 



Where SI*, Sl* A 
the relations 



and S}* x are constants satisfying (E3J) and 



SI* 



1 

37 



SIX 



2 

3t" 



six 



(29) 



Here, E represents a set of infinite saddle points whose 
line element is that of the Einstein universe (k = 1 and 
H = 0). The eigenvalues of these points are determined 
by <j), which in terms of $1* and 7 is given by 



37 
2 



(3 7 - 2) si; + 4 (3 7 - 1) six 



One can check, using J23|) and (psj) , that <j> is always 
positive. The dynamical character of all the equilibrium 
points is given in the table below. As we can see from 
the previous tables, it depends on the equation of state 
(on the parameter 7) and on the expanding (H > <^> 
e = 1) or contracting character (H <0<=>e=— l)of the 



point: 








Model 




Dynamical character 






0< 7 < i 


1=1 


7>| 


F± 


saddle 


saddle 


saddle 


M+ 


repeller 


repeller 


saddle 


M_ 


attractor 


attractor 


saddle 


dS+ 


attractor 


attractor 


attractor 


dS- 


repeller 


repeller 


repeller 


E 




saddle 


saddle 


m+ 


saddle 


repeller 


repeller 


m_ 


saddle 


attractor 


attractor 



At this point, we can observe some differences with the 
general relativistic case ^J. First, the Einstein Universe 
(E) appears to be a critical point for 7 > h, in contrast 
with general relativity, where it appears for 7 > |. On 
the other hand, as we will discuss in detail later, the 
dynamical character of some of the points changes with 
respect to general relativity. For instance, in the brane- 
world scenario the expanding and contracting flat FLRW 
models (F+ and F_ respectively) are no longer repeller 
and attractor, respectively, for 7 > |. They are now 
saddle points. 

Another important difference is that now we have ad- 
ditional critical points, namely, m+ and m_. Let us an- 
alyze in detail the dynamics represented by these mod- 
els. First of all, we have to point out that their char- 
acterization presents an extra difficulty with respect to 
the other models. Their coordinates in the state space 
are SI* = (0,0,0,1) and SI* = (e, 0,0,1), i.e. the con- 
tributions of the ordinary matter term (Sl p ), the spa- 
tial curvature (fife) and the cosmological constant (Sl\) 
are negligible. Therefore, we have at the same time 
n 2 pH~ 2 -> and (6\y 1 n 2 p 2 H^ 2 -> 1, hence their char- 
acterization must involve a limiting process. In order to 
understand the dynamics let us consider the simplified 
situation A = 3 i? = 0, in which the Friedmann equa- 
tion (O) can be solved to give 



tt(t) = (f-M 31 {t + tBB) 3 " 1 , 

where the constant tBB is the Big-Bang time 



(30) 



tBB 



37 2 k 2 A 37 



1_ f« (5) y 
'>7 V K ) 



In the state-space diagrams shown in the Figures be- 
low, this situation corresponds to models in the line join- 
ing m + and F+. From (30), we deduce that for late times, 
t ^S> t BB , the scale factor behaves as a(t) ~ t^y , and 
therefore the solution approaches the flat FLRW model 
(F_|_), hence we have a general relativistic behaviour. 
However, the new interesting behaviour appears when 
we approach the initial singularity (t — + tsB) or, in other 
words, at very high energies (p 3> A), where we have 
a(t) ~ (t — tsB)~ ■ From the point of view of Einstein's 
equations ([|), in such a situation the term involving the 
four-dimensional constant, k , becomes negligible with 
respect to the term involving the five-dimensional one, 
K(5) . We recover general relativity in the limit tBB — > 0, 
which is the opposite situation. From this discussion we 
realize that the limiting process leading to the critical 
points ra± is 



^ 




(31) 



Then, we find that the points m± are models whose scale 
factor is given by 
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l(t) = t* 



(32) 



This is the Binetruy-Deffayet-Langlois (BDL) solution ||] 
(see also |||l0|]) As we have already mentioned, these 
models describe the dynamics near the singularities. 
That is, the early universe behaviour, near the initial 
Big-Bang singularity and also, for recollapsing models 
(for which we must have 3 R > 0), the dynamical be- 
haviour when we approach the Big-Crunch singularity. 
In both cases the dynamics changes with respect to gen- 
eral relativity. 

With the information we have obtained about the crit- 
ical points of the dynamical systems for l~i and H, we can 
apply the well-known techniques used in dynamical sys- 
tems p2| to obtain the structure of the state space, which 
provides, in a visual way, the complete information on 
the evolution of our system (a perfect-fluid FLRW model 
in the brane-world scenario) once the initial conditions 
are given. In the same way as the dynamical character 
of the critical points depend on the equation of state, 
or equivalently, on the parameter 7, so will do the state 
space. In fact, we have found that there are values of 
7 for which bifurcations, that is, topological changes in 
the state space (see |2^] for details), appear. Specifically, 
these values are 73 = 0, |, | (in general relativity we 
only have bifurcations at 7 = 0, |). As we will see in the 
discussion of each particular case, for 7 = -| and 7 = |, 
we have lines with an infinite number of critical points, 
for which we get one vanishing eigenvalue, as is expected 
in those cases [^2| . 

Let us begin with the 7 = case. We have not drawn 
the state space because it is quite simple. Equation (|Io| ) 
implies that the energy density is constant. Then, we can 
solve the Friedmann equation ( |TT| ) and we find that a(t) 
is given by 



a(t) = < 



where U is a constant, e = sgn(iJ), and A is a modified 
cosmological constant given by 



4/4 cosh 






for fc 


= 1, 


< e Vf (*-*«) 




for fc 


= 0, (33) 


A / -2- sinh 
V A 


C\ 


/§(«-**) 


for k 


= -1, 



A = A + ^(l + £ 



(34) 



For A 7^ 0, all the models belong to the de Sitter class, 
whereas in the limit A — > (A = p = 0) we find the 
Minkowski (fc = 0) and Milne (fc = —I) spacetimes. The 
dynamics (of expanding models, e = 1) is reduced to the 
fact that the model fc = is the future attractor, and the 
Milne universe is a repeller. 

For the other cases (7 7^ 0), the whole state space is 
constructed by matching the state space corresponding 
to the dynamical systems (p7-|20) and ( p5f|28| ). It con- 
sists of three pieces, the diagram shown in Figure 



on the right, which corresponds to the case e = 1 in ([l7]- 
^q| ), the diagram in Figure |l](a) in the middle, and on 
the left the diagram corresponding to the case e = — I 
in ( |l7| - |2C| ), which has not been included here because it 
can be obtained from the Figure [l](b) just by reversing 
the direction of the arrows and replacing the subscript 
"+" by "— ". In order to follow the evolution, we have 
specified the quantities represented in the different axes. 
Notice that the state space is compact, with the bound- 
aries given by the planes tt\ = tt\ = 0, Q\ = tl\ = 
and the vacuum models fi p = fl p = 0. 

We have drawn only the trajectories on the planes, 
but the trajectory of any point in the state space outside 
these planes can be deduced qualitatively following the 
behaviour shown in them. As is obvious, the general 
relativistic state space corresponds to the plane fl\ = , 
which is an invariant submanifolc^\ of the state space. 
Therefore, the aim of this work is to study what happens 
when we take initial conditions outside of this plane. The 
other invariant submanifolds are: the vacuum boundary 
fip = 0, the flat geometry submanifold ilk = 0, and the 
SIa = submanifold. 

Keeping this preamble in mind, let us analyze the dif- 
ferent cases according to 7 . For 76 (0, |) and 3 R < 0, 
Milne is a repeller, as in the general-relativistic case, and 
the expanding de Sitter model is the future attractor for 
all the initial conditions excepting the plane = , for 
which the attractor is the flat FLRW model. For 3 R > , 
dS+ plays the same role. In the plane Qa = 0, collapsing 
FLRW models evolve towards the expanding flat FLRW 
model (F + ), with the effect of the extra dimension being 
maximum when H = (<=>• Q = 0). In conclusion, the 
dynamics in this case is essentially the same as in general 
relativity. 

The next case, 7 = g > constitutes a bifurcation. The 
topology of the state space changes [see Figures ||(a) 
and ||(b)] due to the fact that we have now a line of 
vacuum critical points. This line extends to the three 
parts of the whole state space. In the 3 R > sector, 
Figure ||(a), all these critical points, excepting the points 
m± and E, are not included in the previous tables. Their 
coordinates are $7* = (Q*, 0,0,1), \Q*\ < 1, and hence 
they do not appear in general relativity. In order to see to 
what particular models they corresponds we need to con- 
sider the limit ( |3l] ) since they have 0^ = 1. Then, solv- 
ing the Friedmann equation (|ll]), we find they are posi- 
tively curved FLRW models with dynamics described by 
a(t) ~ t. The particular case Q* = corresponds to the 
Einstein universe. In the 3 R < sector, Figure §(b), the 
critical points, excepting points m± and M±, are also not 
in the tables above and they are non-general-relativistic 
in nature. Their coordinates are £1* — (0, 0, ft* x ) with 



* State space trajectories starting in an invariant submanifold 
will never leave it. 
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FIG. 1. State space for the FLRW models with 7 6 (0, ^) and (a) non-negative spatial curvature, 3 R > (on the left) and 
(b) non-positive spatial curvature, 3 R < (on the right). Replacing fl*. by f2 CT and M e by K e the drawing on the right is also 
the state space for Bianchi I models with 7 6 (0, 1) . The planes CIa — Q and S1a — correspond to the state space of general 
relativity. The critical points F e , M e , dS e , E, m e and K e describe the flat FLRW model, the Milne universe, the de Sitter 
model, the Einstein universe, the non-general-relativistic BDL models and the Kasner spacetimes respectively, e is the sign of 
the Hubble function, differentiating between expanding and collapsing models. The planes joining the points dS e , M e and m E 
represent vacuum solutions (fi p = fl p = 0). Only trajectories on the invariant planes, which outline the whole dynamics, are 
drawn (fij. = 0, Oa = £^a = 0, and Q\ = £l\ — 0). 



ill + = 1- Then, using the limiting procedure (31), 
we get the same time dependence: a(t) ~ t . These points 
are also non-general-relativistic. 

The next step is to study the state space in the inter- 
val 7 € (|, |), which is now described by the diagrams 
shown in Figure |. We can find some changes with re- 
spect to the situation in the previous cases. First, we 
have an infinite number of critical points corresponding 
to the Einstein universe, which are arranged in a line de- 
termined by equations (f29|). In the 3 R > sector we 
can see that dS+ and m_ are attractors of the evolution. 
Then, this sector of the state space is divided into two 
regions. The first one consists of those points which will 
evolve to the de Sitter model (dS+), which corresponds 
to the whole 3 R > sector in the case 7 € (0, |). The 
second region is determined by the points which evolve 
towards the BDL model (m_ ) which does not contain any 
general relativistic point. These trajectories correspond 
to models collapsing in the future, that is, evolving to- 
wards a Big Crunch singularity, where the dynamics is 
given by (|32|). It is worth noting that in general rela- 
tivity recollapsing models only occur for 7 > | . These 
two regions are separated by the surface generated by the 
trajectories that start from or arrive to the set of critical 
points representing the Einstein universe (E), which are 
saddle points. In the 3 R < sector the situation is sim- 
pler. For a vanishing cosmological constant = 0) the 
future attractor are the flat FLRW models (F + ) whereas 



in the case with a cosmological constant it is the de Sitter 
model (dS+). 

In 7 = I we have another bifurcation motivated by the 
appearance of two lines of infinite critical points which 
join at the general relativistic Einstein universe, given by 
fl* = (0, 1, 0, 0) (see Figure One of the lines is com- 
posed by Einstein universe points whose state space co- 
ordinates satisfy (p9|). The other line corresponds to gen- 
eral relativistic models (which were not shown in |p7[ ), 
and it occupies the three regions. Their points are char- 
acterized by A = A -1 = , and their scale factor grows 
linearly with time (H ^ 0), a(t) = Ct. They are perfect- 
fluid models with equation of state p + 3p = and energy 
density given by 



2 3(C : 

K p = 



k) 



CH 2 ' 

The case k = — 1 ( 3 R < 0) and C = 1 corresponds to the 
Milne universe. 

The last situation corresponds to the case 7 > |, de- 
scribed by the state space drawn in Figure [| The sit- 
uation in the 3 R > sector is now very similar to that 
showed in the 7 £ (4, | ) case, where two regions ap- 
peared according to whether the points evolve to the 
BDL or to the de Sitter model. The region of points 
evolving to the BDL model is now bigger. With regard 
to the 3 R < sector, the situation has now changed: 
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FIG. 3. State space for the FLRW models with 7 6 (| , | ) and (a) non-negative spatial curvature, 3 R > (on the left) and 
(b) non-positive spatial curvature, 3 7? < (on the right). The drawing on the right is also the state space for Bianchi models 
with 7 G (1, 2). See the caption of Figure |l| for more details. 
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FIG. 4. State space for the FLRW models with 7 = | (a bifurcation) and (a) non-negative spatial curvature, R > (on the 
left) and (b) non-positive spatial curvature, 3 _R < (on the right). The drawing on the right is also the state space for Bianchi 
models with 7 = 2. See the caption of Figure fi] for more details. 



for the models without cosmological constant the attrac- 
tor is the Milne universe (M + ), whereas the flat FLRW 
models (F+) are saddle points. For a non- vanishing cos- 
mological constant, the de Sitter model (dS+) is again 
the attractor. 



III. DYNAMICS OF BIANCHI MODELS IN THE 
BRANE- WORLD SCENARIO 

In this section we will study the dynamics of some ho- 
mogeneous but anisotropic cosmological models (Bianchi 
models) in the brane-world scenario. In particular, we 
will consider the perfect-fluid Bianchi I and V homoge- 
neous cosmological models in which the fluid velocity is 
non-tilted, which means that the hypersurfaces of homo- 
geneity are orthogonal to the fluid flow. Moreover, we 
will also consider a linear equation of state (|J) for the 
perfect fluid. We have considered these two particular 
classes of Bianchi models for simplicity and because they 
contain the flat and negatively curved FLRW models. 

It is well-known that Bianchi models can be described 
by systems of ordinary differential equations, being the 
fluid proper time t the only independent variable that 
appears. The form of the system of ordinary differential 
equations depends on the parametrization of the models, 
i.e., on the variables we use to describe them. Here, we 
will start using the point of view adopted by Ellis and 
MacCallum 129(1 , where they use an orthonormal tetrad, 
{u, e a } (a = 1, . . . , 3), adapted to the fluid velocity 

u ■ u = —1 , u ■ e a = , e a ■ ep = S a p . 



Then, the dynamics can be described in terms of the fol- 
lowing variables: (i) The spatial commutation functions, 
j a ps, defined by the commutation relations between the 
spatial basis vectors, [§p,eg] = ^ a pse a {l a [ps\ = l a ps)- 
Here, we will use the equivalent variables 

u a — 2 ' a P ' — 2 1 ok • 

introduced by Schucking, Kundt and Behr (see |^9| and 
references therein) to decompose "f a ps as follows: "f a ps = 
2a[pS a s] + spg K n aK ■ (h) The kinematical quantities. The 
Hubble function H (= V a u a /3) and the components of 
the shear tensor o~ a b 

°~ab = h a c h b d V \ c u d ) — Hh ab , 

where h a b = 9ab + u a Ub is the orthogonal projector to the 
fluid velocity u. (hi) The matter variables. In our case 
only the energy density p and the isotropic pressure p, 
related by an equation of state (||) . 

In the case of Bianchi models the generalized Fried- 
mann equation reads as follows 

where 2a 2 = a ab a a b and the spatial scalar curvature has 
the following expression in terms of the spatial commu- 
tation functions 

3 R = -6a a a a - n a(3 n a/3 + \ {n a a f . (36) 

On the other hand, from the Einstein equations (||) we 
have a constraint on our variables 



9 




FIG. 5. State space for the FLRW models with 7 > § and (a) non-negative spatial curvature, 3 R > (on the left) and (b) 
non-positive spatial curvature, 3 R < (on the right). See the caption of Figure fyfor more details. 



Zcjupa 13 - e aSK n KT a s T = . 



(37) 



To find systems of equations for Bianchi models similar 
to those described in the FLRW case, we need evolu- 
tion equations for the new variables: a a , n a ® and a a p . 
However, it is better to consider them for each particular 
Bianchi case. 



A. Bianchi I perfect-fluid cosmologies 

The Bianchi I models are homogeneous and anisotropic 
cosmological models containing the flat FLRW space- 
times. We can specialize the triad {e a } in such a way 
that the unit vector fields e a are Fermi- Walker propa- 
gated along u and at the same time their commutation 
functions vanish |29l 



7 5 Q/3 - a a = n a0 = . 



(38) 



Then, in this case the constraint ( |37j ) is identically sat- 
isfied. Moreover, in these models the spatial curvature, 
the curvature of the hypersurfaces orthogonal to the fluid 
velocity, vanishes, that is, 

3 R a b = . 



In particular 3 R = 0, which is a consequence of (36) 
and (|38|). Then, in this case the Friedmann equation (35) 
takes the following form 



n„ + n A + n a + n x = 1 



(39) 



where ft p , fl A and il\ arc defined as in the FLRW case 
[see equations (0,0)], and where we have introduced 



the following dimcnsionlcss quantity associated with the 
shear 



Sin 



3iF 



6H 2 



(40) 



We can construct a state space for the Bianchi I 
cosmological models by taking the variables fi = 
(f2 p , f2A, f2(T, Ha)- Then, taking into account that all 
these quantities are positive by definition, the Friedmann 
equation ( |39| ) implies that we have got a compact state 
space in which these variable are restricted to the inter- 
val [0,1]. Using the time derivative defined in ( |l5| ) and 
using the evolution equation for er 2 p9fl 



(a 2 )' = ~6Ha 2 , 
the system of dynamical equations is given by 
n' p = e[2(l + q)-3 7 }n p , 

n' A = 2e(i + q )n A , 

Q' x = 2e[l + q - 3 7 ]Q X , 



(41) 



(42) 
(43) 
(44) 
(45) 



and the equation (|l6|), which again is uncoupled to the 
rest of equations. Now, the expression for the deceler- 
ation parameter q in terms of the variables H is given 
by 



37- 2 



n A + 2n a + (3 7 - i)n> 



(46) 



The critical points of the dynamical system ([42-45) hav- 



ing a hyperbolic character, together with their state space 
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coordinates f2* = (f2* f2^, f2* , f2^) and their eigenvalues 
are given in the following table: 



dS, 



Model 


Coordinates 


Eigenvalues 




F e 


(1,0,0,0) 


e(3 7 - 2, 37,3(7- 2), - 


-3 7 ) 


dS e 


(0,1,0,0) 


-c(37,2,6,67) 




K e 


(0,0,1,0) 


e(-3( 7 -2),6,4,-6( 7 


-1)) 


m e 


(0,0,0,1) 


e(37,6 7 ,6(7-l), 2(3 7 


-1)) 



where K denotes the Kasner vacuum spacetimes, whose 
line element can be written as follows 



ds^ 



3 

dt 2 + J2 t2pa (dx a f 
=1 



where p a are constants satisfying 



5> Q = - 1 



(47) 



(48) 



a=l 



Apart from the critical points shown in the table above, 
we have found sets of infinite points in the particular 
cases 7 = 0, 1, 2 , which at the same time constitute bi- 
furcations and will be discussed later. From the eigenval- 
ues we get the dynamical character of the critical points, 
which is shown in the table below 



Model 




Dynamical character 






< 7 < 1 


7=1 


7 > 1 


F± 


saddle 


saddle 


saddle 


dS+ 


attractor 


attractor 


attractor 


dS_ 


repeller 


repeller 


repeller 


K+ 


repeller 


repeller 


saddle 


K_ 


attractor 


attractor 


saddle 


m+ 


saddle 


repeller 


repeller 


m_ 


saddle 


attractor 


attractor 



Now let us analyze the state space for these models. 
It can be represented by the same drawings used for the 
3 R < sector of the FLRW evolution, the only thing 
we need to change is the axis corresponding to the vari- 
able 17fc. For the Bianchi I models instead of ilk we have 
to consider CT , and instead of the critical points M± 
(Milne), we have to consider K± (Kasner). Then, tak- 
ing into account this correspondence, let us examine the 
structure of the state space for Bianchi I models for the 
different values of 7. For the sake of simplicity, we will 
do it only for expanding models, that is e = 1. The case 
e = — 1 can be obtaining by a simple time reversal. 

In the case 7 = 0, all the points given by O* = 
(n*,n* A ,0,n* x ) such that Q* + Q* a + fl* x = 1, arc criti- 
cal points corresponding to the de Sitter model as given 
in expressions ([53 34) for k = 0. The dynamics can be 
shown in an one-dimensional state space (see Figure ||) . 
De Sitter is the attractor and Kasner the repeller, there- 
fore, any initial anisotropy is diluted out in the evolution. 

The situation in the case 7 £ (0, 1) is more compli- 
cated. The state space is represented in Figure |l|(b). 
As we can see, F + is the attractor in the case without 
cosmological constant. When a cosmological constant is 



FIG. 6. State space for Bianchi I models with 7 = 0. 



present, the de Sitter is the general attractor. In both 
cases, the anisotropic Kasner models are repellers of the 
evolution, which means that independently of the initial 
conditions the models isotropize. 

When 7 = 1, which corresponds to dust matter, a 
bifurcation occurs and the state space is now given in 
Figure |](b). The dynamical behaviour is essentially the 
same as in the previous case, the difference is that now 
we have an infinite set of critical points situated in the 
line fl a + fl\ = 1 . We can find these models by using the 
limiting procedure introduced in (pl[). The line clement 
of these models is described by the Kasner metric fl47|), 
but now only the first condition in (E8h holds, i.e., 



^2 p ° = 1 ' 



(49) 



where the parameters p a depend on K( 5 ), and in the limit 
p — > 0, where the influence of the extra dimension disap- 
pears, we recover the vacuum Kasner models. 

In the interval 76 (1,2) the only change with respect 
to the situation in the case 7 £ (0, 1) is that now the 
BDL model (m+) is a repeller point and K± are saddle 
points. The state space is represented in Figure ||(b). 
In 7 = 2, which corresponds to a stiff matter equation 
of state (p = p), there is another bifurcation. We have 
a line of general relativistic critical points, as shown by 
Figure |(b). The models are described, as in the pre- 
vious case (7 = 1), by a Kasner metric (E7|) where the 
parameters only satisfy (E^) and the energy density is 
given by 



k 2 p 




(50) 



They are saddle points. The dynamics of the rest of the 
state space is as in the case 7 € (1,2). Finally, for 7 > 2 
(which does not satisfy the causality condition), the state 
space would be given by Figure ^|(b) . 

To sum up, we have seen that expanding models 
isotropize as it happens in general relativity, although 
now we can have intermediate stages in which the 
anisotropy can grow [see Figure |](b) for 7 € (1,2)]. 
The situation near the Big Bang is more interesting. 
In the brane-world scenario anisotropy dominates only 
for 7 < 1, whereas in general relativity it dominates 
for 7 < 2, therefore in the physically relevant interval 
7 G (1,2) the prediction is completely different (for the 
context of inflation see |50| ) : in the brane-world scenario 
the singularity is isotropic. 
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B. Bianchi V perfect-fluid cosmologies 

In the Bianchi V cosmological models the hypersur- 
faces of homogeneity, which we have assumed to be or- 
thogonal to the fluid velocity, are negatively curved. In 
fact, we can pick up a triad {e a } Fermi- Walker prop- 
agated along u and such that the spatial commutation 
functions satisfy 

ai 7^ , a-2 = 03 = , and n a p = . (51) 

Then, the spatial scalar curvature ( |36| ) is given by 

3 R = -6a? < . 



We can introduce the quantity Qk as defined in (12), but 
now it looks as follows 



n, = — 



5 1? 



> 0. 



6H 2 H 2 

Therefore, using also the variables il p , Q\, Q\ [equa- 



tions (y,yj) 

becomes 



and fl„ 



) , the Friedmann equation (|35|) 



Sir 



This equation implies that we can construct a compact 
state space from the variables ft = (Op, £lk, Oa, O ct , fix), 
which are all positive, and as usual, restricted to the 
interval [0, 1]. 

Before looking at the dynamical system for i~2 let us 
consider the constraint (pf\). In this case it is not auto- 
matically satisfied, but it imposes, by virtue of (fill), the 
following condition 



0. 



(52) 



which supposes a restriction on the general metric of the 
Bianchi V models, whose line element can be written in 
the following way 



ds 2 = -dt 2 + A 2 (t)dx 2 + e 2x [B 2 (t)dy 2 
The restriction imposed by (|52|) is then 



C 2 (t)dz 2 



A 1 



BC. 



(53) 



To find the dynamical system we need the evolution 
equations for a\ and a 2 . The equation for a\ is ]29| ] 

di = —Ha\ , 

and the equation for a 2 is ([i"l|). Then, the equations for 
17 are 





= e[2(H 


- q) - 3 7 ]fi p , 


(54) 




= 2egO fe 


5 


(55) 


n' A 


= 2e(l + 


q)Q A , 


(56) 


K 


= 2e(q- 




(57) 


n' x 


= 2e[l + 


q-3j]n x , 


(58) 



where the deceleration parameter is also given by the ex- 
pression (E6|) . The critical points of the dynamical sys- 
tem J54-58) as well as their coordinates and eigenvalues 
are given in the following table |28] 



Model 


Coordinates 


Eigenvalues 


F e 


(1,0,0,0,0) 


e(3 7 - 2, 3 7 -2, 3 7 , 3(7 -2), -3 7 ) 


M e 


(0,1,0,0,0) 


-e(3 7 - 2, 0,-2,4,2(3 7 -l)) 


dS e 


(0,0,1,0,0) 


-e(3 7 ,2,2,6,6 7 ) 


K* 


(0,0,0,1,0) 


e(-3( 7 -2),4,6,4,-6(7-l)) 


m e 


(0,0,0,0,1) 


2e(^, 3 7 -1,3 7 , 3(7-1), 3 7 -l) 



That is, we recover the equilibrium points we had in 
the case of FLRW models with 3 R < plus the mod- 
els denoted by K*, which corresponds to Kasner models. 
However, we must take into account the restriction (|53), 
which implies that the critical points K* only represent 
Kasner models for which the parameters p a are given by 



1 

Pi = g > Vi 



1 + V3 



P3 



1 - \/3 



(59) 



On the other hand, we have now sets of infinite points 
for 7 = 0, 5, |, 1, 2, which also are bifurcation values of 
the parameter 7. We have more bifurcations than in the 
Bianchi I case, so we need more state space diagrams to 
represent the dynamics. To do that we need to extract, 
from the previous table, the dynamical character of the 
equilibrium points, which is shown in the next table 



Model 


< 7 < 1 


Dynamical character 
7 = 1 


7 > 1 


F± 


saddle 


saddle 


saddle 


M± 


saddle 


saddle 


saddle 


dS+ 


attractor 


attractor 


attractor 


dS_ 


repeller 


repeller 


repeller 


K+ 


repeller 


repeller 


saddle 


K* 


attractor 


attractor 


saddle 


m+ 


saddle 


repeller 


repeller 


m_ 


saddle 


attractor 


attractor 



Let us now analyze the state space diagrams shown in 
Figures [^-|lT|. First of all we notice that to get diagrams 
similar to those of the FLRW and Bianchi I models we 
would need four-dimensional diagrams. However, this 
is not necessary since, as it happens in the case of the 
FLRW and Bianchi I models, the qualitative dynamics 
follows from the trajectories of two-dimensional invariant 
submanifolds. Hence, we have drawn two-dimensional 
state space diagrams in which all the dynamical infor- 
mation is present. We have drawn only the trajectories 
joining critical points and the direction of the dynamical 
flow. The interior trajectories can be derived from them 
and by comparison with the state space diagrams for the 
FLRW and Bianchi I models. 

We start with the case 7 = 0, in which the equa- 
tion (|l0|) implies that the energy density is constant. The 
state space is again very simple (see Figure |^), the de Sit- 
ter model, with a modified cosmological constant ([34j), is 
the general attractor and the Kasner and Milne space- 
times are repellers. 
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dS + 



Ml n t 



FIG. 7. State space for Bianchi V models with 7 = 0. 




FIG. 8. State space for Bianchi V models with 7 £ (0, |) . 
Changing the line joining m+ and M+ by a line in which 
all the points are critical we get the state space for 7 = | . 
Reversing the arrow in that line we have the state space for 
7G(§,§). 




FIG. 9. State space for Bianchi V models with 7 = 



a bifurcation characterized by the appearance of a set of 
anisotropic critical points which are not general relativis- 
tic in nature since their coordinates are (0, 0, 0, fi*, 
with f2* + f2^ = 1 . They are among the critical points 
discussed in the case 7 = 1 of Bianchi I. The metric is 
given by the line element (fl7|), but now the exponents p a 
satisfy 



Now, let us consider the case 7 € (0, |) , which is rep- 
resented in Figure ||. Again, for the sake of brevity we 
will consider only the expanding case, i.e., e = 1 (the case 
e = — 1 is obtained by time reversal, i.e. by reversing the 
arrow in the state-space diagrams). The attractor for 
A = is the flat Friedmann model, whereas for A ^ 
is de Sitter. The other critical points are saddle points. 
As in the Bianchi I case, these models isotropize (evolv- 
ing either to F + or dS+), with the exception of extreme 
situations of zero measure, representing models evolving 
towards the BDL solution. For 7 = | we have the first bi- 
furcation due to the appearance of a line with an infinite 
number of critical points located at fl* = (0, fijjl, 0, 0, fij) 



with VL*, 



n* x = 1 



(see Figure H). They are the mod- 



els discussed in the FLRW case (see Figure ||). In the 
case 7 € (|, |) the only change with respect to the case 
7 € (0, ^) is that now models in the line joining m + and 
M+ evolve to M + instead of m + . 

For 7 = 2/3 we have another bifurcation (see Fig- 
ure ^|) with a set of equilibrium points located at 
(fl*, 0, 0, 0), with Q*+Q* k = 1 . They are FLRW mod- 
els that correspond to the critical points in Figure |^(b). 

In the interval 7 G (§,1) (Figure [l0]) we have that 
the flat FLRW model is a saddle point, even if we re- 
strict ourselves to the plane Oa = . For 7 = 1 we have 



1 , 2 

pi = - and p 2 + p 3 = - ■ 



(60) 



In the situation 7 6 (1,2) (Figure [L0|), the only change 
with respect to the case 7 € (§, 1) is that the points 
are saddle instead of repellers. This means that if we 
initially start with f2* + Q* x = 1 , the models will evolve 
towards K^j_ instead of m + . 

Finally, for 7 = 2 we have again an infinite number 
of critical points (see Figure [fl]). They are spatially- 
flat models and therefore, following the discussion of the 
Bianchi I models, they are described by a Kasner met- 
ric ( 47 ) with exponents given by ( |60[ ) and energy density 
by©. 

To summarize, we can say that the dynamics of the 
Bianchi V models encompasses the features of the state 
spaces of the FLRW models with 3 R < and the Bianchi 
I models. 



IV. REMARKS AND CONCLUSIONS 

In this paper we have studied systematically the dy- 
namics of homogeneous cosmological models (the FLRW 
and Bianchi I and V models) in a generalized version 
of the scenario proposed by Randall and Sundrum [0, 
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Kl dS+ K 




FIG. 10. State space for Bianchi V models with 7 6 (|, 1) . 
Changing the lines joining m+ and K!j_ by lines in which all the 
points are critical we get the state space for 7 = 1. Reversing 
the arrow in those lines we have the state space for 7 G (1, 2) . 



K* dS+ k; 




FIG. 11. State space for Bianchi V models with 7 = 2. 



stressing the main differences with respect to the general- 
relativistic case. In the case of the FLRW cosmological 
models, the state space presents a new equilibrium point, 
namely, the BDL model (m±) ||. It dominates the dy- 
namics at high energies, where the extra-dimension ef- 
fects become dominant. For this reason, we expect them 
to be a generic feature of the state space of more gen- 
eral cosmological models in the brane-world scenario, as 
it occurs in the Bianchi models analyzed here. In the 
FLRW case the critical points m_|- describe the new dy- 
namics near the Big Bang and also near the Big Crunch 
for recollapsing models. 

Another new feature is the existence of new bifurca- 
tions as we change the equation of state, the parameter 
7 . In the case of FLRW models there is one new bifur- 
cation for 7 — i , characterized by the appearance of an 
infinite number of non-general-relativistic critical points. 
Among them we find a static model whose line-element is 
that of the Einstein universe. This contrasts with general 
relativity, where it appears for 7 = | . The consequence 
is that in the brane-world scenario recollapsing models 
appear for 7 > -| instead of 7 > | as in general rela- 
tivity. In the case of Bianchi I models we have found a 
new bifurcation for 7 = 1, and in the case of Bianchi V 
models for 7 = i , 1 . 

On the other hand, Bianchi models allow us to study 
anisotropy. We have seen that expanding Bianchi I and 
V models always isotropize, as it happens in general rel- 
ativity, although now we can have intermediate stages 
in which the anisotropy grows. This is an expected re- 
sult since the energy density decreases and hence, the 
effect of the extra dimension becomes less and less im- 
portant. The situation changes drastically when we look 
backwards. Near the Big Bang the anisotropy only domi- 
nates for 7 < 1, whereas in general relativity it dominates 
for 7 < 2, which includes all the physically interesting 
cases. 

Just to finish we would like to mention some current 
and future work in the line of the present one. First we 
recall that in this paper we have considered brane-world 
scenarios in which the bulk satisfies the condition (Q). 
Then, it would be interesting to look at the effect of hav- 
ing a contribution from the bulk curvature, or in other 
words, a contribution from the bulk Weyl tensor piece 
E^ h . This is currently under investigation |2(J. On 
the other hand, taking into account that string theory 
is formulated in spacetimes with more than one extra di- 
mension (brane world of codimcnsion greater than one) 
it would be interesting to study how the introduction 
of more extra dimensions changes the results presented 
here. In this sense, a good starting point would be to 
consider scenarios like those introduced in |B1|. 
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